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EQUILIBRIUM GEOCHEMICAL MODELING OF A SEASONAL THERMAL ENERGY 

STORAGE AQUIFER FIELD TEST 

J. A. Stottlemyre 
Pacific Northwest Laboratory 


ABSTRACT 


The report sunmarizes a geochemical mathematical modeling study designed 
to investigate the well plugging problems encountered at the Auburn University 
experimental field tests. The results, primarily of qualitative interest, 
include: 1) loss of injectivity was probably due to a combination of native 
particulate plugging and clay swelling and dispersion, 2) fluid-fluid 
incompatabilities, hydrothermal reactions, and oxidation reactions were of 
insignificant magnitude or too slow to have contributed markedly to the 
plugging, and 3) the potential for and contributions from temperature-induced 
dissolved gas solubility reductions, capillary boundary layer viscosity 
increases, and microstructural deformation cannot be deconvolved from the 
available data. 


INTRODUCTION 


Since 1976, two field experiments have been conducted by Auburn Univer- 
sity involving storage of heated waters in a shallow, confined aquifer near 
Mobile, Alabama (ref. 1 and 2). The first experiment involved injection of 
7,570 m^ of 37°C, filtered water from an electric power plant cooling water 
canal. The storage aquifer is located between 40 and 62 m below the land 
surface. The heated canal water was stored for 37 days and then recovered 
with an overall thermal efficiency of 67%. The process was inhibited by 
significant plugging of the injection well. This has been attributed to clay 
and silt particles suspended in the canal water; filtering above the 5-micron 
range improved but did not eliminate the problem. 

A second experiment, utilizing the same storage aquifer, involved a 79-day 
injection of 55,345 m^ of 55°C boiler heated water. The water source was an 
unconfined aquifer located between 25 m and 34 m below the land surface. The 
water was stored for 50 days and then retrieved over a 41-day period. The 
recovery efficiency was 65% over a temperature range of 55°C and 3°C. Ambient 
groundwater temperature was 20°C. Figure 1 is a schematic of this second 
experiment. It is important to note that the supply water was extracted from 
an overlying aquifer and, therefore, the system did not represent a true 
doublet configuration. 

Clogging of the injection well again proved to be a major operational 
difficulty. Loss of permeability resulted in a decrease in the maximum 
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injection rate from 12.6 f.»sec"l (200 gpm) to 6.3 £«sec“l (100 gpm) . 

Plugging of the well may have been due to the water sensitive nature of the 
storage aquifer sediments. Mon tmoril Ionite clay in combination with low 
cation concentration of the supply water relative to the storage aquifer water 
may have resulted in swelling and dispersion of clay particles as shown con- 
ceptually in figure 2. Such water sensitivity is a documented phenomenon which 
lends itself to laboratory identification and field pretreatment (ref. 3, 4, 

5). It is also possible that the supply aquifer water contained suspended 
solids and/or dissolved gases which may have contributed to the plugging. 

In general, there are other potential reservoir permeability damage mecha- 
nisms including precipitation of minerals due to the mixing of incompatible 
groundwaters, water-rock incompatibility, increased temperatures, boundary 
layer viscosity anomalies (ref. 6), and microstructural deformation (ref. 7). 

The primary objective of this study is to investigate fluid-fluid incom- 
patibility, fluid-rock incompatibility, hydrothermal mineral alterations, and 
redox reactions with respect to potential contribution to the loss of well 
injectivity observed at the Auburn field experiments. This investigation was 
based on equilibrium chemical thermodynamic computer modeling. No laboratory 
and/or post -experimental field data are available for comparison, and there- 
fore the results of this computer study are only of qualitative value. 


Sediment and Groundwater Characterization 

Approximate groundwater chemistry is shown in table 1. Sediment min- 
eralogy and grain size distribution are given in tables 2 and 3. 


Description of Test Cases 

The objective of this study is to analyze some potential alternate causes 
of formation plugging at the Auburn field test. The following observations 
are noted : 


• The shallow partially confined supply aquifer is low in ionic concentra- 
tion relative to the storage aquifer, and has unknown suspended solid and 
dissolved gas concentrations. 

• The deeper, confined, storage aquifer water is assumed to be in chemical 
equilibrium with the formation mineralogy at a temperature of 19.5°C. 

• The supply and storage waters were increased in temperature from 19. 5° C 
to approximately 55°C. 
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• Plugging occurred when the supply aquifer water was used as the working 
fluid. 

• Plugging apparently did not occur when the storage formation water was 
injected back into the storage aquifer. This water was not heated however 
(Molz, Auburn University, personal communication). 

Based on these observations and a water sensitivity test conducted at 
Auburn University, (Molz, personal communication), it is probable that the loss 
in well injectivity resulted primarily from particle plugging and clay swelling 
and dispersion. The montmoril Ionite content is less than a percent by weight; 
however, this is often sufficient to inhibit the flow of low salinity waters 
(ref. 3). Furthermore, dispersed particles in the micron and submicron range 
can often significantly reduce permeability (ref. 5). 

Other potential reservoir damage mechanisms amenable to chemical 
thermodynamic modeling include: 

• mineral precipitation as the working fluid temperature is increased 

• mineral precipitation as oxygen is introduced to the system 

• mineral precipitation due to the mixing of the potentially incompatible 
supply and storage formation waters 

• mineral precipitation due to a chemical incompatibility between the sup- 
ply water and the storage aquifer sediment and/or hydrothermal alteration 
products . 

To study these four potential categories, several computer simulations 
were conducted as listed in table 4. 


Data Analysis and Conclusions 


For each computer simulation, the following equilibrium data were tabu- 
lated: 1) equilibrium mineralogy, 2) type and quantity of new mineral precipi- 

tates, 3) fluid temperature and pH, and 4) aqueous species concentrations. It 
is assumed that if minerals precipitate to any significant degree, decreased 
formation permeability might result. 

Increasing the working fluid, temperature and/or oxygen content (Eh) appar- 
ently has a negligible effect on mineral precipitation. Hematite, a ferric 
(iron) oxide, is the only mineral susceptible to precipitation. However, as 
shown in table 5, the quantity in moles per kilogram of water is rather insig- 
nificant. In addition, mixing the supply aquifer water and the storage aquifer 
water does not result in deleterious mineral precipitation; therefore, fluid- 
fluid incompatibility should be discounted as a contributing factor in the 
observed formation plugging at the Auburn field tests. 
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Alteration of the storage aquifer mineralogy was also investigated. This 
scenario involves interacting the heated supply water with ambient temperature 
groundwater and sediments in the storage aquifer. Table 6 shows the ionic con- 
centrations in the groundwater and equilibrium mineralogy predicted by EQUILIB 
for four different injection water temperatures. With respect to the min- 
eralogy at 55°C, it can be observed that EQUILIB predicts that calcite, musco- 
vite, and kaolinite all react to some extent and that the minerals adularia 
and calcium montmorillonite would be formed as products. Similarly, the feld- 
spars mi croc line and low-albite apparently react and zoisite is a predicted 
reaction product. As the injection water temperature is increased, there is a 
net decrease in the amount of solid material within any volume of rock equili- 
brated at these temperatures. To maintain a mass balance, there is an increase 
in the aqueous species concentration of the fluid. This might indicate that 
as the temperature decreases with increasing distance from the well, precipi- 
tation may occur. If the fluid is saturated with respect to certain mineral 
species at elevated temperature near the well, transport of the fluid to a 
lower temperature environment could result in precipitate formation. The con- 
sequences of such precipitation would depend, in part, on the quantity and 
density of the precipitate and the interstitial makeup of the sedimentary 
matrix. However, kinetics is an additional factor that must be considered. 

Based on equilibrium predictions, it might be argued that hydrothermal 
mineral alteration contributed to the plugging observed at the Auburn test 
site. However, the computer results should be viewed with caution and consid- 
ered to be qualitative only. Because of thermodynamic inconsistencies in the 
data base and the equilibrium assumption, results predicted by complex geo- 
chemical computer codes may not always be accurate. The assumption that the 
heated supply water, the storage aquifer water, and the storage aquifer min- 
eralogy have achieved a stable equilibrium becomes quite restrictive if the 
temperature under consideration is as low as 55°C. Reaction rates of rock- 
forming minerals with aqueous solutions may be extremely slow. Hydrothermal 
reactions probably do not occur rapidly enough to account for the plugging 
observed in the first 48 hours of the Auburn experiment. 

It is concluded that with the possible exception of clay swelling and 
dispersion, fluid-rock incompatibility was not a contributing factor in the 
formation damage in the Auburn tests. Furthermore, heating the supply water, 
introducing oxygen, and mixing the supply and storage formation waters appar- 
ently had no effect on precipitation of minerals or the creation of alteration 
products that could reasonably explain the formation damage. Based on the 
limited evidence, it is assumed that water sensitivity (clay swelling and dis- 
persion), particulate plugging, and outgassing of dissolved gasses represent 
the most reasonable explanations for the loss of injectivity in this specific 
case. More detailed study might reveal that temperature-induced capillary 
boundary layer viscosity anomalies or microstructural deformation may have been 
contributing factors. Table 7 is a summary of the potential formation damage 
mechanisms and a qualitative estimate of the relative likelihood of each having 
occurred at the Auburn facilities. 
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TABLE 1. GROUNDWATER CHEMISTRIES AS SUPPLIED BY AUBURN UNIVERSITY 


Supply Aquifer Storage Aquifer 

(Unconfined Aquifer) (Confined Aquifer) 

mg/& (mg/& 


Water Type 

Sample 1 

Sample 2 

Sample 1 

Sample 2 

Na 

3.40 

0.60 

7.60 

11.10 

Ca 

0.33 

0.21 

0.05 

0.38 

Fe 

0.07 

0.07 

0.06 

0.05 

Si 

6.30 

9.10 

10.40 

9.70 

CaCO 

3.80 

3.80 

176.00 

176.00 

pH 

7.19 

7.19 

7.38 

7.38 

Temperature 

19.5°C 

19.5°C 

19.5°C 

19.5°C 


TABLE 2. INITIAL SAMPLE MINERALOGY BASED ON OPTICAL PETROGRAPHY 
AND X-RAY DIFFRACTION 

Mineral Composition Concentration (Vol%) 


Calcite 

CaC0 3 

3.7 

Quartz 

Si0 2 

76.5 

Hematite 

Fe 2°3 

2.5 

Muscovite 

KAl 2 (AlSi 3 )0 1Q (0H) 2 

1.3 

Kaolinite 

Al 4 Si 4°10 (OH) 8 

2.8 

Alkali Feldspar 

(K,Na)AlSi 3 0 g 

8.3 

Plagioclase Feldspar 

NaAlSi 3 0g - CaAl 2 Si 2 0g 

4.5 

Montmoril Ionite 

(l/2Ca,Na,K) 0 7 (Al,Mg,Fe) 4 



(Al,Si) g O 20 (OH) 4 • nH 2 0 

0.4 


614 



TABLE 3. GRAIN SIZE DISTRIBUTION FROM SIEVE ANALYSIS 


Size Fraction Height Percent 


18 x 

35 

0.22 

35 x 

120 

86.39 

120 x 

200 

4.82 

200 x 

325 

4.18 

-325 


4.39 


Description 

Phi Size 

Coarse Sand 

0. 0-1.0 

Medium Sand 
Fine Sand 

1. 0-3.0 

Very Fine Sand 

3.0-3. 7 

Very Fine Sand 
Coarse Silt 

3. 7-4.5 

Coarse Silt 
and Finer 



TABLE 4. COMPUTER TEST CASES 

Input Oxygen 20 °C 55°C 100°C 150°C 

1) Supply Water Alone - - X X X 

2) Supply Hater Alone X X X X X 

3) Storage Water Alone - X X X X 

4) Supply Water Plus Storage 

Water X X X X X 


5) Supply Water Plus Storage 

Water Plus Minerals X - X X X 
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TABLE 5. MINERALS PRECIPITATED DUE TO HEAT, OXYGEN AND/OR FLUID-FLUID INCOMPATIBILITY 


Fluid 

Oxygen 

Temperature (°C) 

pH 

Eh 

Insoluble Minerals 

Quantity 
(moles/kg water) 

Supply Water 
Alone 

Yes 

20 

7.19 

+0.807 

quartz Si 02 
hematite Fe 203 

0. 192E-03 
0.627E-06 



55 

8.36 

+0.541 

hematite Fe 203 

0.627E-06 . 



100 

7.77 

+0.354 

hematite Fe 203 

0.627E-06 



150 

7.34 

+0.226 

hematite Fe 203 

0.627E-06 

Supply Water 
Alone 

No 

55 

8.36 

-0.215 

minnesotaite Fe 3 Si 40 io(OH )2 

0.418E-06 ; 



100 

7.77 

-0.159 

hematite Fe 203 

0.627E-06 



150 

7.34 

-0.178 

hematite Fe 203 

0.627E-06 

Storage Water 
Alone 

No 

20 

7.38 

-0.091 

quartz Si 02 

minnesotaite Fe 3 Si 40 i(j( 0 H )2 

0.274E-03) 

0.328E-06 



55 

8.95 

-0.268 

minnesotaite Fe 3 Si 40 io(OH )2 

0.328E-06 



100 

8.36 

-0.196 

hematite Fe 203 

0.492E-06 



150 

7.90 

-0.212 

hematite Fe 203 

0.492E-06 

Supply Water 
Plus Storage 
Water 

Yes 

20 

7.27 

+0.367 

quartz Si 02 
hematite Fe 203 

0.234E-03 

0.560E-06 



55 

8.76 

-0.079 

hematite Fe 203 

0.560E-06 



100 

8.17 

-0.120 

hematite Fe 203 

0.560E-06 



150 

7.71 

-0.151 

hematite Fe 203 

0.560E-06 



TABLE 6. WATER AND MINERAL EQUILIBRIUM CONCENTRATIONS AT VARIOUS TEMPERATURES 


Moles/kg H7O 


Aqueous Species 

20 °C 

55°C 

100°C 

150°C 

N + 

0.247-03 

0.125-02 

0.125-02 

0.125-02 

Ca ++ 

0.611-05 

0.614-04 

0.284-03 

0.601-03 

Fe ++ 

0.492-06 

0.318-13 

0.193-22 

0.288-16 

Fe +++ 

0.627-06 

0.405-13 

0.270-12 

0.547-11 

Si 

0.316-03 

0.291-03 

0.845-03 

0.215-02 

CO3 - 

0.132-03 

0.476-02 

0.407-02 

0.828-02 

0H“ 

0.100-06 

- 

- 

- 

H + 

0.100-05 

- 

- 

- 

0 2 

0.390-03 

- 

- 

- 

Ai+++ 

- 

0.124-05 

0.203-05 

0.720-05 

K+ 

- 

0.690-03 

0.234-02 

0.233-02 

pH 

7.27 

8.260 

7.270 

6.760 

Eh 

- 

+0.018 

+0.769 

+0.361 

p C0 

- 

- 

0.047 

0.417 

Mineralogy 

Calcite (CaC03) 

0.100-01 

0.537-02 

0.606-02 

0.186-02 

Quartz (Si0 2 ) 

0.400-02 

0.609-02 

0.505-02 

0.751-02 

Hematite (Fe 2 03> 

0.100-02 

0.100-02 

0.100-02 

0.100-02 

Muscovite (KAl3Si302o(OH) 2 ) 

0.100-02 

0.280-05 

- 

- 

Kaolinite (Al^Si^OjoCOHig) 

0.100-02 

0.707-03 

- 

— 

Microcline (KAlSigOg) 

0.100-02 

- 

- 

- 

Low Albite (NaAlSi30g) 

k-Montmorillonite Kq. y(Al,Mg,Fe)4(AI,Si)g0 2 o(OH)4 

0.100-02 


" 

" 

• nH 2 0 

0.100-02 

- 

- 

— 

Adularia (KAlSigOg) 

- 

0.164-02 

- 

- 

Zoisite (Ca 2 Al3Si 3 0 12 (0H)) 

Ca-Montmorillonite (l/2Ca)o. 7(Al,Mg,Fe)4(Al, Si)g 


0.228-02 

0.160-02 

0.377-02 

0 2 q(OH )4 • nH 2 0 

- 

- 

0.280-02 

- 

Total 

2.00-02 

1.710-02 

1.65-02 

1.41-02 



TABLE 7. POTENTIAL DAMAGE MECHANISMS FOR THE AUBURN FIELD TESTS 


Mechanism 


Qualitative 

Potential Comments 


Temperature- Induced Phenomena 

• Mineral precipitation Low 

• Out gas sing Moderate 


• Increased quartz-water Unknown 

boundary layer viscosity 

• Micros true tural deformation Unknown 


Fluid-Fluid and Fluid-Rock 


Chemical Incompatibility 

• Clay swelling and High 

dispersion 

• Mineral precipitation Low 

(fluid mixing) 

• Mineral precipitation Low 

(oxidation) 


Fluid-Rock Physical Incompatibility 
• Suspended solids Moderate 


• Existing formation 


High 


Mathematical modeling potential 
for dissolved oxygen in the 
supply aquifer 
Limited available data 

Limited available data 


Significant montmorillonite, 
low salinity water injection 

Mathematical modeling 

Mathematical modeling 


Potential for micron and sub- 
micron particles in the supply 
aquifer water 

Loose clay and silt particles 
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FIGURE 1. SCHEMATIC DIAGRAM OF THE SECOND MOBILE, ALABAMA, 
FIELD TEST SYSTEM (ref. 2) 
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FIGURE 2. CONCEPTUALIZATION OF CLAY SWELLING AND DISPERSION 
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